Size-dependent melting: Numerical calculations of the phonon spectrum 
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In order to clarify the relationship between the phonon spectra of nanoparticles and their melting 
temperature, we studied in detail the size-dependent low energy vibration modes. A minimum 
model with atoms on a lattice and harmonic potentials for neighboring atoms is used to reveal a 
general behavior. By calculating the phonon spectra for a series of nanoparticles of two lattice types 
in different sizes, we found that density of low energy modes increases as the size of nanoparticles 
decreases, and this density increasing causes decreasing of melting temperature. Size-dependent 
behavior of the phonon spectra accounts for typical properties of surface-premelting and irregular 
melting temperature on fine scales. These results show that our minimum model captures main 
physics of nanoparticles. Therefore, more physical characteristics for nanoparticles of certain types 
can be given by phonons and microscopic potential models. 



I. INTRODUCTION 

The physical properties of nanosized material have at- 
tracted considerable attention because of its technologi- 
cal importance as well as its fundamental interest for the- 
oretical study. One important thermodynamic property 
is its size-dependent melting point depression, which, in 
fact, has been studied for decades. This unusual behav- 
ior of nanoparticles can lead to develop new features of 
materials with wide applicationsii^i^. 

Experimental studies show that (with a few excep- 
tions), for free surface nanoparticles, the melting tem- 
perature decreases with decreasing size, and has irregu- 
lar variations on a fine scale^i^. Many theoretical mod- 
els have been employed to investigate this unusual prop- 
erty. Of these various models we can classify them into 
many categories, thermodynamical models^i^iiS.iiiilS,, la- 
tent heat modeli^ii^, the surface-phonon instability 
modeli^ii^i^, the liquid drop modeli^i-, and bond en- 
ergy model22ii^. Most of these models give the result 
that, the nanoparticle melting temperature Tmn relative 
to the bulk melting temperature Tmb decreases linearly 
with size in a general form such as Tmn/Tmb = 1 — -^1 D, 
where D is the nanoparticle diameter and A is a constant 
depending upon the modeling parameters. Some other 
models have a nonlinear form to fit better for smaller 
sizes {D < 10nm)^i22. 

The melting properties of two-dimensional Penrose lat- 
tices have been studied by Huang et. al, who made use 
of the phonon spectrum and the Lindemann melting 
criterion^^. It showed that the boundary atoms always 
have larger atomic mean-squared thermal vibration am- 
plitudes, which means lower melting temperature. We 
will see that a similar behavior also appears in three- 
dimensional lattices in our calculation. 

In this paper, we study the relation between nanoparti- 



cle melting temperature and nanoparticle size by means 
of the phonon spectrum, and we intend to reveal the 
melting properties of nanoparticles from a microscopic 
model. With a statistical melting criterion similar to the 
Lindemann criterion, we found that the linear depression 
of nanoparticle melting temperature with decreasing size 
and the irregular variations on a fine scale was given in 
our numerical calculation by a very simple model. 



II. MODEL 

Two three-dimensional lattice types are considered 
in this paper. For simplicity, we only calculate near- 
spherical nanoparticles in this paper. To construct such 
nanoparticles, first, we create the simple cubic (SC) or 
face-centered cubic (FCC) lattices with n lattice points 
on the edge. Then cut an inscribed sphere in the corre- 
sponding lattice type. The lattice points on and inside 
the sphere are what we count in the calculation (actually 
there are no lattice points on the sphere for even n SC 
lattice). Assume that there is one atom at each lattice 
point, which is the equilibrium position of the atom, and 
the surface is free. So we do not have to take surface 
reconstruction and relaxation into account in this calcu- 
lation. 

We directly calculate the phonon modes in the har- 
monic approximation for the two lattices of different sizes 
(different n). The interatomic potentials that we use is 
the spring potential 
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where is the position of the ith lattice point and it' is 
the displacement of the atom from its equilibrium posi- 
tion Tq-, (ij) means that the sum is taken over the near- 
est neighbors and {ij} means the next-nearest neighbors; 
Ki and K2 are spring constants for the nearest and next- 
nearest interactions respectively. This potential ensures 
that the structure of nanoparticle we consider is at least 
one of the lowest energy structures. 

Within the harmonic approximation, we can deduce 
that the atomic mean-squared thermal vibration ampli- 
tude (u2^) is 

{^l) = j:^i^jn^ + h (2) 

k 

where i — 1,2, ... ,N and N is the number of atoms 
of the nanoparticle; a = x,y, z; rrii is the mass of the 
ith atom; lj^ is the fcth eigenfrequency, and Cia^k is the 
component of the eigenvector corresponding to the fcth 
eigenfrequency. We set h = ks = 1 here. ( ) means grand 
canonical ensemble average. 

In bulk materials, the Lindemann criterion^^ is com- 
monly used. It states that a material melts at the tem- 
perature at which the amplitude of thermal vibration ex- 
ceeds a certain critical fraction of the interatomic dis- 
tance. Experiments and simulations show that its criti- 
cal value is around 0.10-0.15 in units of the interatomic 
distance. Lawson showed that the overall accuracy of 
the Lindemann criterion for the materials he considered 
is not better than about 30%, and he gave an improved 
Lindemann criterion^^. For irregular nanoparticles, the 
distance-fluctuation criterion is introduced^^i^. It is 
based on the fluctuation of the distance between pairs 
of atoms, while the Lindemann criterion is based on the 
fluctuation of individual atoms relative to their average 
position. And it has a problem that it depends on the 
duration of the simulation ru n^^i^^ . 

It can be seen from the expression of (uf^) that, when 
the atoms all have the same mass rrii = m, the to- 
tal mean-squared atomic thermal vibration amplitude 
{u'^)tot = ^^("D is determined only by the eigenfre- 

quencies except the mass m, because of the orthonor- 
mality of the eigenvectors ^ s^^^k = 1' i-^- 
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In this case, we can directly use the Lindemann criterion 
in our model and do not need to calculate the eigenvec- 
tors at all. 

But when the atoms do not have the same mass, we 
have to and actually we can easily calculate each {u^^) 
numerically here. Therefor it is reasonable to use a statis- 
tical criterion like the Lindemann criterion: iJ more than 
fifty percent of the \/ (uf^) exceed half of the interatomic 
distance, the melting happens. The "fifty" and "half are 



not particularly selected and can be tuned. Compared 
to the Lindemann criterion, the effect of using these two 
figures together may be too large. But obviously they 
can be tuned to be consistent with the Lindemann cri- 
terion. And because the model we use here is so sim- 
ple and we just focus on qualitative description, there 
should be no problem to use these two figures in the mo- 
ment. We will show that, with this statistical criterion, 
the linear depression of nanoparticle melting tempera- 
ture with decreasing size and the irregular variations on 
a fine scale are reproduced in our model and a crosspoint 
in the phonon modes and related results are discovered 
when we consider {uf^) statistically, although we only 
consider the case that all the atoms have the same mass. 

Actually, the eigenfrequencies depend only on the mass 
rrii and the spring constant Ki, K2 in our model. So each 
at a given temperature is completely determined 
when rrii and Ki , K2 have been chosen. Different lattice 
spacing a has different melting temperature because 
of its appearance as equilibrium interatomic distance in 
the criterion. 



III. RESULTS AND DISCUSSION 

In the present calculation, we only consider the case 
of monoatomic nanoparticles. The SC type lattice has 
the nearest-neighbor interaction as well as the next- 
nearest-neighbor interactions, while FCC type only has 
the nearest-neighbor interaction. We set m.i = m ~ 
1, Ki = 3, K2 = 2, a = 1 for SC and to^ = to = 1, Ki = 
9, i^z = 0, a = for FCC. We choose a = 1 for SC 
and a — -s/A for FCC so that they have the same number 
density. 

The density of phonon modes for lattices in two dif- 
ferent finite sizes along with their thermodynamic limit 
are shown in Fig. [1] for SC and FCC respectively. When 
the size of nanoparticle becomes larger, the profile of the 
density of phonon modes becomes smoother. 

For the two finite sizes of SC lattices, there are both 
three peaks at w « 2.4, 3.4, and 4.6. The peak at w « 
2.4 lowers down with increasing size, and it becomes a 
shoulder at a higher frequency lo w 2.8 ultimately in the 
thermodynamic limit. So does the peak at w « 4.6, but 
it increases with increasing size and the shoulder in the 
thermodynamic limit is at a lower frequency uj « 4.4. 
The peak aX lo k, 3.4 for finite sizes becomes the only 
peak at a higher frequency lo ~ 3.7 in the thermodynamic 
limit. This peak in the thermodynamic limit arises from 
the symmetry that leads to a large number of degenerate 
eigenfrequencies at this value. And our calculation show 
that when the symmetry breaks, the degeneracy is lifted 
and peaks beside it appear in finite sizes. The situation 
for FCC lattices is similar. But there are no shoulders in 
the thermodynamic limit and the peak at a; w 3.5 in the 
smaller finite size disappear in the larger finite size while 
the peak at w w 7.5 does not change much for FCC. 

The density of phonon modes for nanoparticles of those 
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FIG. 1: (Color online) Density of phonon modes for different 
sizes, (a) SC lattices. The dashed (blue) line is for n = 13, 
the dotted (red) line is for n — 20, and the solid (green) line 
is for n — > oo. (b) FCC lattices. The dashed (blue) line is 
for n = 9, the dotted (red) line is for n = 13, and the solid 
(green) line is for n —f oo. n is the number of lattice points on 
the edge of the cubic lattice, from which we cut the inscribed 
sphere. 

sizes, which are smaller than the smaller finite size (n — 
13 for SC and n — % for FCC) appears in Fig.[l] begins to 
have zigzag-like profiles and to fiuctuate more and more 
strongly when the size gets smaller and smaller. 

To see which phonon modes contribute to the melt- 
ing, we pick out those {uf^) that meet our criterion and 
plot the mean contribution from each [w, w-j-dw] to {uf^) 
versus uj at the nanoparticle melting temperature deter- 
mined by our model. The mean contribution from each 
[u, Lo + Aijj] to (m^q.) is defined as 

where means the sum is taken over those {uf^) that 
meet our criterion, and 

^ — ' rriiLOk e^''! ^ — 1 I 

k 

means the sum of k in the expression of (uf^) (Eq. ^ 
is taken only over those k that are in the range of ujk G 
[liJ,(jJ + dw], and Nm means the number of those (uf^). 
We choose duj = 0.2 for SC and dcu = 0.4 for FCC. 

Four different finite sizes for SC and FCC lattices are 
shown in Fig. [5] respectively. Noted that there is a cross 
point in both parts of the figure, i.e. uj « 2.6 for SC and 
uj « 4.5 for FCC, so the eigenfrequencies are divided into 
two parts naturally: low and high eigenfrequencies. For 
the smallest size (the dotted lines) shown in Fig. O most 
contributions are from the low eigenfrequencies and the 
line fluctuates strongly in the low eigenfrequency range. 
With the increase of the size, the contributions from the 
high eigenfrequencies increase, but the increase becomes 
smaller and smaller; the contributions from the low eigen- 
frequencies get smaller with increasing size and also there 
is a smaller and smaller decrease, but the change is much 
bigger than that of the high eigenfrequencies. Appar- 
ently, a peak corresponding to the peak at w w 2.4 for SC 
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FIG. 2: (Color online) Mean contributions from different 
eigenfrequencies to the melting, (a) SC lattices. The dot- 
ted (brown) line is for n = 5, the dotdashed (green) line is for 
n = 9, the dashed (blue) line is for n = 13, and the solid (red) 
line is for n — 20. (b) FCC lattices. The dotted (brown) line 
is are for n = 4, the dotdashed (green) line is for n = 6, the 
dashed (blue) line is for n = 8, and the solid (red) line is for 
n = 13. n is the number of lattice points on the edge of the 
cubic lattice, from which we cut the inscribed sphere. 

and at « 3.5 for FCC in Fig. [T] respectively appears in 
the lines (the dotdashed line and the dashed line) for the 
two middle sizes. Similarly, a peak corresponding to the 
peak at w « 4.6 for SC and at w ~ 7.5 for FCC appears 
in the line (the solid line) for the largest size in Fig. ^ 
which will be apparent if we plot it alone. So the line of 
contribution becomes more similar to the line of density 
of phonon modes, when the size becomes larger. The con- 
tributions almost vanish at the highest eigenfrequencies. 
So the low eigenfrequencies dominate the melting for rel- 
atively small nanoparticles and they can be excited at low 
temperature. But for relatively large nanoparticles, the 
high eigenfrequencies nearly give the same contributions 
to the melting. That is why the nanoparticle melting 
temperature decreases with decreasing size. 

We can show that all the phonon modes with small 
mean neighbor relative deviations, ((u* — u^Y), belong 
to the low eigenfrequency part. We call these acoustic- 
like phonon modes. Most phonon modes with large mean 
neighbor relative deviations are within the high eigen- 
frequency part, and we call them optical-like phonon 
modes. In the acoustic-like modes, neighbor atoms tend 
to vibrate in the same direction similar to the acoustic 
modes in the diatomic linear chain; while in the optical- 
like modes, they prefer to vibrate in the opposite direc- 
tion. So with the same thermal vibration energy, the 
acoustic-like modes can result in larger atomic mean- 
squared thermal vibration amplitude, but the optical-like 
modes can not do that because they have large potential 
energy when the neighbor atoms move to each other with 
small displacements. 

Then we divide the nanoparticle into five layers for SC 
and four layers for FCC, each layer having a mean ra- 
dius r from the center of the inscribed sphere, and plot 
the contributions from the low and high eigenfrequen- 
cies respectively versus r (i.e. P{r) = -^^py^ Z]"("L)Aa;, 
where means the sum is taken over those i, a whose 
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corresponding atoms are in the same layer and {uf^)Auj 
means the sum of k in the expression of (uf^) is taken 
over those k whose corresponding ujk belong to the low or 
high eigenfrequencies; (u'^)tot appears for normalizaion) 
in Fig.[3]for SC and Fig.Hfor FCC, both with four differ- 
ent finite sizes respectively. Here all {uf^) are considered 
no matter whether they contribute to the melting or not. 
This is different from Fig. [21 

Apparently, the low eigenfrequencies vibrate the atoms 
near the surface most, and they have nearly no effect on 
the inner atoms. It means surface-premelting, because 
the low eigenfrequencies are excited at low temperature 
and they mainly vibrate the atoms near the surface. This 
behavior should be a general property in the model simi- 
lar to ours. The high eigenfrequencies vibrate the atoms 
in the middle when the size is relatively small, and they 
vibrate the atoms all but the inner ones when the size be- 
comes large. So when most atoms are near the surface, 
the low eigenfrequencies is sufficient for melting because 
of their low excitation temperature, and when the size 
becomes large, more and more atoms are in the middle, 
the high eigenfrequencies have to be excited at high tem- 
perature for melting. 

The relation between nanoparticle melting tempera- 
ture and nanoparticle size are shown in Fig. [5l The 
data calculated for SC and FCC lattices are fitted by 
Tmb — A' /D. We do not use I — A/D because we have no 
idea on the bulk melting temperature Tmb for our model. 
It can be seen that the linear decrease of the nanoparticle 
melting temperature with decreasing size is reproduced 
for the two lattices. And the nanoparticle melting tem- 
perature fluctuates stronger with decreasing size. We 
think the reason is that when the size decreases, the den- 
sity of phonon modes gets irregular, so it is uneasy to 
find a regular line to fit these small sizes. A = A' /Tmb 
is approximately 2.0 for SC and 2.5 for FCC. The dif- 
ference comes from the fluctuation of the nanoparticle 
melting temperature of small sizes. If we only consider 
the large sizes calculated in our model, SC and FCC will 
have nearly the same A in our calculation {A = 2.2 for 
SC and A = 2.3 for FCC). 



IV. CONCLUSIONS 

In this paper, we numerically calculate the phonon 
spectra for two types of three-dimensional lattices (SC 
and FCC), in the framework of a very simple microscopic 
model with a spring interatomic interaction within the 
harmonic approximation. Each spatial component of the 
mean-squared thermal vibration amplitude of each atom 
can be obtained during the calculation, and we are able 
to introduce a statistical melting criterion similar to the 
Lindemann criterion. Within this criterion, the linear 
decrease of the nanoparticle melting temperature with 
decreasing nanoparticle size and irregular variations on 
a fine scale are reproduced by using our simple model. 
And we found that the eigenfrequencies are naturally di- 
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FIG. 3: (Color online) Contributions from the low and high 
eigenfrequencies from five layers with mean radius r from the 
center of the inscribed sphere for SC lattices. The circles 
(red) are for low eigenfrequencies and the squares (blue) are 
for high eigenfrequencies. (a) n = 5 (b) n = 9 (c) n = 13 (d) 
n = 20. n is the number of lattice points on the edge of the 
cubic lattice, from which we cut the inscribed sphere. 
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FIG. 4; (Color online) Contributions from the low and high 
eigenfrequencies from four layers with mean radius r from the 
center of the inscribed sphere for FCC lattices. The circles 
(red) are for the low eigenfrequencies and the squares (blue) 
are for the high eigenfrequencies. (a) n = 4 (b) n = 6 (c) 
n — 9 (d) n = 13. n is the number of lattice points on the 
edge of the cubic lattice, from which we cut the inscribed 
sphere. 
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FIG. 5: (Color online)The relation between nanoparticle 
melting temperature and nanoparticle size, (a) SC lattice, 
n is from 5 to 20. The fit line is T = 2.03 - 4.06536/D. 
(b) FCC lattice, n is from 4 to 13. The fit line is T = 
2.05495 - 5.07076/i3. 



vided into two parts, the low eigenfrequency part and 



high eigenfrequency part, each of which ought to in- 
clude all the acoustic-like and most of the optical-like 
modes respectively. The low eigenfrequency part played 
a major role in the melting of small nanoparticles, which 
resulted in the depression of the melting temperature. 
The atoms mostly near the surface are mainly vibrated 
in these modes, that is the so-called surface-premelting. 
This should be a general feature when considering the 
mean-squared thermal vibration amplitude of each atom 
independently. In this paper, we only consider the case 
of monoatomic nanoparticles, and the case of diatomic 
nanoparticles is under investigation. Moreover, we are 
planning to employ a more realistic interatomic interac- 
tion (including the electron-phonon coupling) and intro- 
duce relaxation of the structure in the further study. 
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